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ABSTRACT 

In the inner regions of an accretion disk around a black hole, relativistic protons 
can interact with ambient matter to produce electrons, positrons and 7-rays. The re- 
sultant steady state electron and positron particle distributions are self-consistently 
computed taking into account Coulomb and Compton cooling, e^e"*" pair production 
(due to 7 — 7 annihilation) and pair annihilation. While earlier works used the diffu- 
sion approximation to obtain the particle distributions, here we solve a more general 
integro-differential equation that correctly takes into account the large change in par- 
ticle energy that occur when the leptons Compton scatter off hard X-rays. Thus this 
formalism can also be applied to the hard state of black hole systems, where the dom- 
inant ambient photons are hard X-rays. The corresponding photon energy spectrum 
is calculated and compared with broadband data of black hole binaries in different 
spectral states. The results indicate that the 7-ray spectra {E > 0.8 MeV) of both the 
soft and hard spectral states and the entire hard X-ray/7-ray spectrum of the ultra- 
soft state, could be due to p — p interactions. These results are consistent with the 
hypothesis that there always exists in these systems a 7-ray spectral component due 
to p — p interactions which can contribute between 0.5 to 10% of the total bolometric 
luminosty. The model predicts that GLAST would be able to detect black hole bina- 
ries and provide evidence for the presence of non-thermal protons which in turn would 
give insight into the energy dissipation process and jet formation in these systems. 

Key vifords: accretion, accretion disks — black hole physics 



1 INTRODUCTION 

Black hole X-ray binaries are generally observed to be in 
two distinct states. These states, which are named hard and 
soft, differ in their luminosity and spectral shapes. In the 
hard state, which is in general less luminous, the spectrum 
of the system can be described as a hard power-law with 
a spectral index F ~ 1.7 and a cutoff around 100 keV. 
This spectrum can be modeled as thermal Comptonization 
of soft p hotons by a plasma h aving temperature T ~ 50 
keV (e.g. iGierlinski et alJllQQTi) . In contrast, the spectrum 
during the soft state consists of a blackbody-like component 
with kT ~ 1 keV, which typically dominates the luminos- 
ity and is generally considered to be thermal emission from 
an optically thick accretion disk. Apart from this soft (or 
disk) component, a hard X-ray power-law tail, with photon 
index F ■-^ 2.5 and no detectable cutoff up to ~ 8 MeV 
iMcConnell et aDl2002^ . is also observed. 
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There have been several interpretations of the high 
energy {E > 200 keV) emission from black hole sys- 
tems. It may arise from a photon starved inner most re- 
gion of a disk which cools due to bremsstrahlung self- 
Comptonization (iMelia Misralll993f) or as emission due 
to tt'' decay, which are created by pr oton-proton interaction 
in a hot prot on gas, T > 10^^ K iKolvkhalov fc SunvaevI 
I1979I: Jourdai nfc Roauejll99l . While these possibilities 
maybe viable, they are based on assumptions of the ge- 
ometry and physical properties of the system, which are 
not directly and independently verifiable, like the presence 
of a very hot proton gas or a photon starved region. An- 
other interpretation is that the spectra arises due to Comp- 
tonization of photon s by the bulk motion of ma tter falling 
into the black hole jLau rent fc Titarchu3ll999l) . However, 
iNiedzwiecki fc Zdziarskin200^ have argued that the non- 
detection of spectral breaks at _E < 500 keV is contrary to 
this model's prediction. A detailed radiative model, which 
has been used to fit good quality broad band data of black 
hole systems, is the hybrid model which is inscribed in 
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a spectral fittinR code call ed EQPAIR jPoutanen fc Coppil 
ll998l:lGierliriski et alJl999l) . In the framework of this model, 
this high energy component arises from a plasma consisting 
of both thermal and non-thermal electrons that Comptonize 
external soft photons. Detailed spectral modeling of the ob- 
served soft state spectra demands the coexistence of thermal 
and non-thermal electrons in the emission region and hence 
the steady state non-thermal electron distribution is com- 
puted by assuming that there is an injection of non-thermal 
particles into the system where they cool by Comptonization 
and Coulomb interactions. 

The origin of these non-thermal particles is uncertain. 
A possible site may be a corona on top o f a cold disk which is 
heated by magnetic field reconn ections jHaardt fc Maraschil 
Il993l: iPoutanen fc FabianlllQQSl) . but the details of the pro- 
cess are largely unknown. The acceleration process has 
to be highly efficient to produce non-thermal electrons 
in an environment where electrons cool rapidly by in- 
verse Comptonization. If this acceleration process is mass- 
independent then protons are also expected to be ac- 
celerated to relativistic energies, for example by scatter- 
ing off magnetic "kinks" in a Keple rian accretion disk 
iSubramanian. Becker fc Kafatos^ll996^ . Some of these non- 
thermal protons may escape from the system a nd contribute 
to the jet formation iSubramanian et alJll999l) . This is par- 
ticularly interesting since there are some evidence that the 
X-ray pro ducing region may b e same as the base of the ex- 
tended jet iMarkoff et al.l2005^ . Hence the detection of these 
non-thermal protons will provide valuable clues to the na- 
ture of black hole systems. 

Non-thermal protons would interact with the ambient 
thermal protons and produce electron-positron pairs, which 
would Comptonize photons to high energies. These high en- 
ergy photons would produce further pairs by 7 — 7 inter- 
action and a pair cascade would ensue. Pair cascades ini- 
tiated by the injection of pairs (or equivalently high en- 
ergy non-thermal electrons) have been extensively examined 
(e.g. iLightman fc Zdziarskil [19871 ISvenssonI [l987l) . The ef- 



fect oi p — p interactions and the resul tant spectra have 



also been computed and studied 
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eral, did not consider the presence of copious photons or 
have assumed that the ambient photons are in the UV 
range, a scenario relevant to AGN and under-luminous black 
hole systems. However, for black hole binaries, the system 
is dominated by either soft or hard X-rays depending on 
the spectral state. One of the important radiative interac- 
tion that would occur in such an environment is the inverse 
Compton scattering of X-rays by pairs with Lorentz fac- 
tor 7 ~ 200 that are produced by the p ~ p interaction. 
The standard methods to compute the inverse Compton 
spectra, assume that the interaction in the rest frame of 
the electron takes place in the non-relativistic Thompson 
limit, which is true only if electron Lorentz factor times the 
photon energy 76 <C meC^. This assumption is violated if 
the ambient photon energy is ^ rrieC^ /"f « 2 keV. Thus, 
when Bhattacharvva ct al. (2003) considered p — p interac- 
tion in the presence of blackbody photons having tempera- 
ture T « 1 keV, they used a general formalism to describe 
the in verse Compton process given bv lBlumenthal fc Gouldl 
il970h . They found that for such a situation, which is rele- 



vant to the soft state of black hole binaries, the effect of non- 
thermal protons is to produce a broad feature around 1 — 50 
MeV. Using the observed OSSE data for GRS 1915-^105, 
they could constrain the fraction of non-thermal protons 
in the system to be < 5%. Although iBhattacharvva et alJ 
i2003h . computed the change in energy of the photon (and 
hence the change in energy of the lepton) appropriately in 
the Klein-Nishina regime, they used, for simplicity, a diffu- 
sion equation to describe the kinetic evolution of the pairs, 
which intrinsically assumes that the change in energy of 
the particle per scattering is small. This assumption breaks 
down when the ambient photon energy is in X-rays, espe- 
cially when there are a copious amount of ambient hard 
X-rays. 

In this work, we e xtend the formalism developed 
by IBhattacharvva et al.l i2003l) . by solving an integro- 
differential equation for the pair kinetic evolution which cor- 
rectly takes into account the large energy changes that a lep- 
ton undergoes upon scattering with an X-ray photon. This 
not only allows for a more accurate estimation of the emer- 
gent spectra for systems in the soft state, but enables the 
scheme to be applied to the hard state also. 

In §2 we describe the model and the assumptions made 
to compute the steady state non-thermal electron/positron 
distributions and the resultant photon spectra. In §3 some 
general results of the computation are presented along with 
comparison with observations of black hole systems Cyg X- 
1 and GRS 1915-f 105 in different spectral states. The main 
results of the work are summarized and discussed in §4. 



2 PAIR DISTRIBUTION AND RADIATION 
SPECTRA COMPUTATION 

We consider a uniform sphere of non-relativistic thermal 
plasma with number density nr and radius R, in the pres- 
ence of an external copious photon source. It is convenient to 
parameterize the luminosity of the external photons Lph, in 
terms of the compactness parameter lph = Lpha-T/{RmeC^), 
where ctt is the Thomson cross-section. The spectral shape 
of the ambient photons is taken either to be a Wien peak or 
an exponentially cutoff power-law depending on the spectral 
state being considered. In this system, we assume that there 
is a power-law distribution of non-thermal protons with in- 
dex a, which would lead to proton-proton collisions i.e. the 
density of non-thermal protons is given by 

nNT{'y)d'y = no{-y - iy°'dj (1) 

The normalization Uo of this distribution is characterized by 
the compactness parameter Ip-p = Lp-par / {RnieC?), where 
Lp-p is the total power in electron, positrons and 7-rays 
which would be produced in the proton-proton interactions. 

The steady state positron density N+ (7) is determined 
by solving the integro-differential equation 

|-(7ciV+(7)) + iV+(7) d7'P(7,7') - 

poo 

/ d7'P(7',7)iV+(7')+iV+(7) = Q+,pp(7)+Q+.^^(7) (2) 
J ~i 
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while the corresponding electron density N-{-y) is obtained 
from 

|-(7civ_(7)) + iv_(7) y"\yp(7,y)- 



d7'P(7',7)iV_(7') =0-.pp(7) + Q-,^^(7) 



(3) 



Here, -P(7, 7 — e)dedt is the probability that a 
positron/electron with Lorentz factor 7 will suffer a colli- 
sion and its Lorentz factor changes by an amount between 
e and e + de in time dt, where e = hv/m^c? is the nor- 
malized photon energy. Q±,-y7 and Q±,pp are the creation 
rates of pairs due to p — p interactions and photon-photon 
production respectively, N+{'y) is the annihilation rates of 
positrons with the thermal background electrons and 7c 
is the Coulomb cooling rate. This formalism to obtain the 
particle distribution follows from the results described in 
iBlumenthal &: Gouldl il970l) (and references therein), where 
the exact expression for P(7,7 — e)dedt has been derived. 
For small changes in particle en ergy, eqs. and can be 
reduced to diffusion equations jBlumenthal fc Gouldlll97dl 
i.e. the integral terms can be reduced to ^ ['y/c-^i (7)] 
where 7/0 is the inverse Compton cool ing rate . It was these 
diffusi on equations that were used bv iBhattacharvva et alJ 
l)20n.^ as a si mplifying assumpt ion. 

Following lEilek fc Kafatod lll983h . e^e~ production 
rate due to p — p process is given by 



:(7p)/3pW]vt(7p) 



(7* 



mil' 



I 3/2^ 

+ 7pfc ; 



1/2 



where pp = (1 — l/7p)^''^ and ^pk = 1 — 7p. The denominator 
of the integrand represents the appropriate energy distribu- 
tion of the leptons produced, while the limits impose the 
allowed energy range and are given by 

JvAl) = 1 + [7*7 + (7' - If^'il' - 1)'/' - 1]'^' (5) 



and 



7p,i(7) = 1 + [7*7 - (jt - 1)^''^(7^ - 1) 



\l/2 



1] 



4/3 



(6) 



where 7* = 70. The approximate but analytical cross- 
section (a^±) is tabul ated for different energy ranges by 
lEilek fc Kafato^ ^8^. 

Pair production rate from photon-photon interactions 
is approximated to be 



Q±,i-i{'y) = c / ".7(27e — e)n^(e)cr^^(e, 27 - e)de 



(7) 



where it has been assumed that for two photons annihi- 
lating with energies e and e', the resultant Lorentz factor 
is ~ (e + e')/2. The app roximate fo rm for the cross sec- 
tion a^f^(e,e ) is given bv ICoppi fclB landford (199(|. The 
positrons primarily annihilate with the background thermal 
electrons at a rate given by 



^+(7) = N+{'y)nTa^+^- {1,^)0 



(8) 



where the approximate for m of the cross section is given by 
ICoppi fc BlandfordI lll99d) . 

It is worthwhile to note that, while the numerical com- 
putation of the integro-differential equations (eqs.|5|and0| 
are relatively straight forward, there are a few minor but 



tricky points that have to be taken into consideration. In 
particular, the number conservation of photons for the in- 
verse Compton process has to be strictly satisfied in the 
numerical computation, in order to obtain the correct and 
stable particle density solution. 

The equilibrium photon density inside the sphere is a 
solution of 



(9) 



r c 

L H[l + rKN(e)\ 

where R-^~f{e) is the rate of photon annihilation and tkn is 
the Klein-Nishina optical depth. Q-^jc, Qj.pp and Q^,e+e- 
are the photon production rates due to inverse Compton, p — 
p interaction and pair annihilation respectively. The photon 
production due to p ~ p reaction is given by 



(^TT- i'yp)l3pnNT {lp)d^p 



(10) 



For simplicity, the scattering of high energy photons 
with the thermal particles are neglected here and the impli- 
cations of this assumption are discussed in the last section. 
The interactions of photons with high energy protons, which 
requires a much larger threshold prot on energy (7^ > 300), 
have also been neglected. As noted bv IBhattacharvva et alJ 
i2003h . such interactions are not important when a > 2 and 
depend on the unknown high energy cutoff of the proton 
acceleration process. 

Equations Q and @ are solved self-consistently 
to obtain electron and positron distributions as well as the 
radiative flux. The output spectrum depends on the spectral 
shape of the external photons and three other parameters: 
the Thomson optical depth r, photon compactness Iph and 
the ratio of proton- proton compactness Ip-p, to that of the 
external photon, (3 = Ip-p/lph- The results are insensitive 
to the non-thermal proton index a and apart for an overall 
normalization factor, to the size R of the system. From the 
definition of compactness {L oc IR) , it follows that for a fixed 
value of I, the normalization of the photon spectrum will be 
proportional to R. The spectral shape primarily depends on 
the inverse Compton and pair production processes, which 
are characterized by tic ~ n~f^s(JTR and r^^ ~ u^cttR- 
Since the soft photon density n^,^ oc Lph/R? oc Iph/R and 
similarly, the high energy photon density, cx Ip-p/R, tic 
and T-y^ are independent of R for specified compactness. 
Thus, parameterizing the luminosities in terms of compact- 
ness renders the shape of the spectrum to be nearly inde- 
pendent of the size of the system. 



3 RESULTS 
3.1 Soft State 

The computed radiated spectra due to p — p interactions, 
corresponding to parameters typical of the soft state of black 
hole binaries, are shown in Figure Since the motivation 
of this work is to investigate the possibility that the 7-ray 
part of the spectrum is due to p — p interactions, a detailed 
fit to the broad band data, especially the low energy part of 
the spectrum, is not being attempted. Instead, the dominant 
soft component is being represented by a a Wien peak like 
spectrum at kTph = 1 keV and high compactness, Ipu = 250. 
The corresponding luminosity is 
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Figure 1. The radiated spectra for parameters typical of the 
soft state of black hole binaries, which are Iph = 250, kTph = 1 
keV, T = 2.5 and a = 2.5. The three curves labeled as 1, 2 
and 3 are for compactness ratio /3 = Ip—p/lph = 0.001, 0.01 and 
0.1 respectively. The solid lines represent spectra obtained using 
the integro-differential approach to solve for the particle energy 
distribution. The dotted lines represent spectra obtained using 
the diffusion approximation as described in FBhattacharvva ct al. 




(^) (_JL.) ,„) 

where R ~ 10^ cm is the size of the system. The optical 
depth of the thermal electrons is of order unity and is taken 
to be T = 2.5 in Figure Q The spectra plotted in the fig- 
ure, correspond to three different values of the compactness 
ratio /3 = Ip-p/lph ~ 0.001, 0.01 and 0.1, which represent 
increasing efficiency of the acceleration process to produce 
non-thermal protons. At low /3 values the spectra are hard 
with a break around ~ 80 MeV. This cutoff is due to the 
interaction of higher energy photons with the ambient soft 
photons of energy ~ 3 keV to generate pairs. The hard spec- 
tra are due to inverse Comptonization of photons by a steady 
state pair distribution N±('y) cc 7"^ (Figure Such a dis- 
tribution is expected when pairs are injected at a large en- 
ergy (7 ~ 200) and are cooled by inverse Comptonization. 
Essentially in the diffusion approximation, the cooling rate 
71C oc 7^, which leads to A'^(7) oc 7"^ when pair creation due 
to photon-photon interaction can be neglected. For higher 
compactness ratio, /3 = 0.01, pair creation due to photon- 
photon interaction affects the particle distribution at low 7 
(Figurel^J, which leads to a softer output spectra (Figure^. 
For still larger values of /3 = 0.1, apart from the change in 
particle distribution due to photon-photon pair production, 
another break at around 511 keV appears in the radiated 
spectrum, since the density of high energy photons is large 
enough for pairs to be produced by self interaction. 

Figure0also shows for comparison, the computed spec- 
tra obtained using the diffusion approximatio n to solve for 
the pa rticle distribution as described in Bhatta charvva et alJ 
i2003r . While, quantitative differences can be seen when the 
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Figure 2. The particle distributions corresponding to the spec- 
tra shown in Figure Q. The dotted (solid) lines represent the 
positron (electron) distributions. 



compactness ratio is large (spectra marked 2 and 3), the 
qualitative features of the spectra are similar. 

The high energy (5-50) soft state spectra are charac- 
terized by a steep power-law of photon index F ~ 2.5, and 
hence cannot be explained as emission from p — p interac- 
tions, since in this energy range the expected spectral index 
is F < 2 (Figure^. Indeed, the soft state spectra can be de- 
scribed by the hybrid model, where non-thermal electrons 
are injected into a plasma. In Figurej^we reproduce the un- 
folded data points based on such a model (the "EQPAIR" 
model) as fitted bv lMcConnell et al.l (|20o3). These data ob- 
tained from instruments on board the BeppoSAX and CGRO 
satellites cover a wide energ y range from ~ 0.5 keV to ~ 8 
MeV. As shown bv .McConneU et al.l (|20o3), the spectrum 
at energies < 1 MeV can be represented as an extension 
of the non-thermal spectra at hard X-rays. On the other 
hand, the 7-ray (i.e. the 1-8 MeV) spectrum, maybe due 
to an additional component, especially since there are un- 
certainties about the relative normalization of the different 
instruments used to make the composite spectrum. A second 
component interpretation is supported by the non-detection 
by COMPTEL of any flux in the energy range from 750 keV 
to 1 MeV, indicating a possible spectral break in that re- 
gion, although the low statistics of the data does not allow 
for a definite conclusion. Thus, we attempt to fit the 7-ray 
part of the spectrum, as emission from p — p interactions, 
us ing para meters th at are consistent with those obtained 
bv lMcConnell et alJ ll2002l) . The required compactness ratio 
P = Ip-p/lph ~ 0.046, can be compared with Int/lph = 0.11 
obtained by the EQPAIR fit. This indicates that if the 7- 
ray spectrum is due to p — p collisions, the energetics of 
the proton acceleration that produces non-thermal protons 
is comparable to the electron one. The factor of two dif- 
ference, between the two ratios, could be either due to the 
neutrino emission, which has not been added to Ip-p or to 
the uncertainties involved in the theoretical modeling and 



Self- consistent computation of'y-ray spectra due to proton-proton interactions in black hole systems 




E(keV) 



Figure 3. Broad band s oft state data of Cyg X-1 reproduced from 
iMcConnell et aTl (2223) • Solid line is the computed spectrum for 
/3 = 0.046, Iph = 125, kTph = 0.4 keV, t = 2.5 and a = 2.5. 
The thick broad hne shows the estimated GLAST sensitivity 
: / /www-glast. sIac.stanford.edu/software/IS / glast_lati_perfoi 



fitting of the data. If the 7 ray spectrum is not due to this 
process, the above estimate serves as an upper limit on the 
energetics of the non-thermal acceleration of the protons. 
For the parameters used to compute the p-p component, 
the ratio of the non-thermal proton density to the thermal 
one turns out to be ~ 0.03, which is consistent with that 
esti mated earlier by Bhatta charyya ct al. ( 200^). 

IZdziarski fc Gierliiiskil i2004l) have estimated the sensi- 
tivity of a 4 X 10^ s GLAST observation to be « 10"^ keV 
cm~^ s~^ at 10 GeV. Figure|3]shows that the predicted spec- 
trum at that energy is a couple of order of magnitude higher 
than this limit and would be easily detecta ble. The extension 
of the hybrid plasma model fit to the data ijMcConnell et all 
|2002^ to 10 GeV also pred icts a similar flux. However, as 
pointed out by |zdziarski fc Gierliriskil i2004) . the hybrid 
plasma model assumes a low value of compactness l^h ~ 4 
and a high value of the maximum Lorentz factor of the elec- 
trons, 7mai = 10*. For a more realistic value of compactness 
or if electrons are not accelerated to such high energies, there 
will be a sharp cutoff in the hybrid plasma spectrum at « 1 
GeV. Thus, a GLAST observation of Cyg X-1 (or black hole 
binaries in general) during the soft state, should be able to 
discern between these models. 

Black hole systems exhibit a variation of the soft state 
spectrum, which is sometimes denoted as a separate state 
and is called Very High or Intermediate state (VHS/IS). 
Here the hard X-ray power-law is nearly equally luminous as 
the soft component. The black hole system GRS 19154-105 
displays a wide variety of spectral shape and its x class be- 
havior (as classified by jBelloni et. aL ( 200(1) l is similar to 
the V HS state of other black hole systems. IZdziarski et alJ 
l|200lfl have fitted the broad band data of this source ob- 
tained from instruments on board RXTE and GGRO, by the 
hybrid model and the unfolded data based on that model is 
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Figure 4. Broad band Ver y-High state (y class) data of GRS 
1915-fl05 reproduced from IZdziarski et alj 1200 if) . Solid line is 
computed spectrum for f3 = 0.075, Ipfi = 50, kTpi^ = 1.25 keV, 
T = 2. 5 and a = 2.5 
e.htm| . 

reproduced in Figure |1| The absence of high energy COMP- 
TEL data does not allow for detailed fits, but a comparison 
of the 600 keV flux with the computed spectrum provides 
an upper limit of the compactness ratio /? = 0.075, which 
again can be compared to the ratio Int/hh = 0.11 required 
by the hybrid model fit JZdziarski et alJl200ll) . 

Another variation of the soft state spectra is the so 
called ultra-soft state, where the spectrum is dominated by 
the soft blackbody like emission, with a weak hard X-ray tail. 
This state has been observed in many black hole systems like 
GRS 19154-105, GX 339-4 and XTE J1550-564. An example 
of this stat e is shown in F i gure El Here the data is repro- 
duced from IZdziarski et al.l i200ll) . and is unfolded based on 
the EQPAIR model fit to observation of GRS 1915-1-105 in 
the 7 class using RXTE and CGRO satellites. A remark- 
able observational feature is that for all these sources, the 
weak high energy tail is hard with a photon spectral index 
r ~ 2, for which currently there is no theoretical explana- 
tion iZdziarski fc Gierliiiski2004ll . For high compactness the 
spectra due to p — p interaction, is expected to have such a 
hard spectral slope (Figure 0. Indeed as shown in Figure |H| 
the hard X-ray spectra can be explained solely as emission 
from such a process with a compactness ratio j3 ~ 0.008. 
Again, the absence of data greater than 1 MeV does not al- 
low for a clear identification. However, the data is consistent 
with the hypothesis, that in the Ultra-soft state, the emis- 
sion from accelerated non-thermal electrons is absent and 
hence, only the hard spectra arising from p ~ p interactions 
is observed. 



3.2 Hard State 

Figures [H] and Q show the computed spectra and particle 
densities for different values of compactness ratio /3, when 
the ambient photon density is similar to that found in the 
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Figure 5. Broad band Ult ra-soft state (7 class) data of GRS 
1915+105 reproduced from IZdziarski et alj fcoOl) Solid line is 
computed spectrum for /3 = 0.008, Iph = 555, kTpi^ = 1.35 keV, 
T = 5.0 and a = 2.5 
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Figure 7. The particle distributions corresponding to the spec- 
tra shown in Figure j^. The dotted (solid) lines represent the 
positron (electron) distributions. 
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Figure 6. The radiated spectra for parameters typical of the hard 
state of black hole binaries, which are Ipf^ = 250, Ec = 100 keV, 
T = 1.0 and a = 2.5. The three curves labeled as 1, 2 and 3 are 
for compactness ratio /3 = Ip—p/lph = 0.001, 0.01 and 0.1 respec- 
tively. The solid lines represent the total spectra, the dotted lines 
represent the spectral component due to p— p interactions and the 
dot dashed line represents the thermal component. The dashed 
line represent the spectral compon ents obtained using the diff u- 
sion approximation as described in lBhattacharvva et al.| ^2003^ 



hard state of black hole binaries. In this state, the spectrum 
is dominated by a thermal Comptonized component, which 
is approximated here as an exponentially cut off power law 
with photon index T = 1.7 and cut off energy Ec = 100 
keV. As the average photon energy in the plasma is higher 
than that of the soft state, the cut-off in the photon spec- 
trum due to pair production appears at « 3 MeV instead 
of « 80 MeV. In the soft state, the shape of the computed 
spectra and particle density depended on the compactness 
ratio p. This was due to increasing self-interaction of the 
non-thermal photons to produce pairs as compactness in- 
creased. However, in the hard state, since the number den- 
sity of ambient photons with energy > 500 keV is large, 
pair production is dominated by interaction of non-thermal 
photons with the background ambient ones even for large 
compactness. Thus, the shape of the computed spectra is 
relatively invariant to l3 as shown in Figure |S] and in partic- 
ular the slope of the expected spectrum at energies > MeV 
is insensitive to the parameters of the model. 

Figure|S|also shows for comparison, the computed spec- 
tral components due to p — p interaction, obtained using the 
diffusion approxim ation to solve for the particle distribu- 
tion as described in iBhattacharvva" et al.l lipoi). Similar to 
what was found for the soft state (Figure 0, only quanti- 
tative differences can be seen while the qualitative features 
of the spectra are similar. This suggests that even when the 
ambient photons have energy ~ 100 keV, the diffusion ap- 
proximation can reveal the salient features of the high energy 
spectrum. 

The broadband hard state data of Cyg X-1, using 
RXTE and BeppoSAX satellites, were analyzed and fitted 
using the EQPAIR model by McConncU ct al. ( 2002). The 
unfolded data is reproduced in Figure|H| Most of the spectra 
can be represented by thermal Comptonization, but there is 
evidence for a different high energy spectral component in 
the 0.7—8 MeV range, which was explained as emission from 
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Figu re 8. Broad band hard state data of Cyg X-1 reproduced 
from iMcConnell et al.l (2223) ■ Solid line is computed spectrum 
for /3 = 0.018, Iph = 110, Ec = 120 keV, r = 1.0 and a = 2.5. 
The thick broad line shows the estimated GLAST sensitivity 
jhttp://www-glast. slac.stanford.edu/software/IS/glast_lati-perfor 
The dotted line represents the spectral component due to p — jj 
interactions. 

non-thermal electrons by iMcConnell et all l|20o3) (see also 
llbragimov et al.l (l2005ll ). Here we compare the data with the 
computed spectra due to p — p interactions, using parame- 
ters representative of the hard state and show that the the 
MeV data, especially the spectral slope in that energy range, 
match with the observations. The required compactness ra- 
tio P = 0.018 is similar to the values required to fit the soft 
state data, b ut smaller than Int/lh = 0.08 required by the 
EQPAIR fit jMcConnell et alJl2002l) . The overall computed 
spectra fits even the data at energies < 700 keV remarkably 
well, however, this could be a coincidence given that an ex- 
ponentially cutoff power-law is not a good approximation 
to a Comptonized thermal spectrum and that reflection is 
not taken into account here. Nevertheless, the fit in Figure 
|5] strongly indicates that the emission from Cyg X-1 in the 
hard state can be explained as thermal Comptonization and 
a component due to p — p interactions. 

The spectrum shown in Figure |S1 predicts that GLAST 
should be able to detect Cyg X-1 in the hard state also. 
Similar to the soft state, if the high energy emission is due 
to non-thermal e l ectron s (like in the EQPAIR model fit by 
IMcConnell et al.l (|20o3)) a detection by GLAST would only 
be possible under extreme conditions like low compactness 
and h igh maximum Lorentz factor of the non-thermal elec- 
trons iZdziarski fc Gierliriskl2004l) . Thus, Gi^ST would be 
able to differentiate between these models. 



4 SUMMARY AND CONCLUSION 

A self-consistent scheme, which computes the elec- 
tron/positron and radiation energy distribution inside a 
thermal plasma having non-thermal protons, is developed. 



The non-thermal protons interact via p — p collisions to pro- 
duce electron/positron pairs and 7-rays. These high energy 
pairs cool by inverse Comptonization of ambient photons. 
The effect of subsequent pair production due to photon- 
photon interactions and annihilation of positrons with ambi- 
ent electrons, is taken into account. Unlike previous works, 
the scheme does not assume that the energy change of the 
leptons during inverse Compton scattering is small. Hence 
it can be applied to black hole binary systems, where the 
dominant ambient photons are X-rays. 

Comparison of the computed spectra with the broad 
band X-ray/7-ray spectra of black hole binaries in different 
spectral states reveal: 

• The hard X-ray spectra (3-800 keV) of the soft and Very 
High state (VHS), are too steep to be explained as emission 
due to p — p interactions. However, the observed 7-ray spec- 
tra (0.8-8 MeV), especially for the soft state, could be due 
to such interactions. In this interpretation, an unknown ac- 
celeration process energises both electron and protons, pro- 
ducing their non-thermal distributions. While the electron 
non-thermal distribution gives rise to the hard X-ray emis- 
sion, the protons are the origin of the 7-rays. The powers 
going into accelerating protons and electrons are within a 
factor of two. 

Vnr the ultra-soft state, the observed hard X-ray spectra 
^.h^L J- j^gy-j gQyijj ]-,g entirely due to p — p interactions. 
This interpretation gives a natural explanation for the sim- 
ilar hard X-ray spectral slope (F ~ 2) observed whenever a 
black hole is in the ultra-soft state. 

• For the hard state, the observed 7-ray spectrum (0.5-8 
MeV) can be explained as emission due to p — p interac- 
tions alone. The predicted steep spectral shape in this en- 
ergy range, is not sensitive to the model parameters, and 
matches well with the observations. 

• For both the soft and hard states, the model predicts that 
for reasonable parameters, GLAST should be able to detect 
black hole binaries. This in contrast to the situation when 
only non-thermal electrons are present in the system, where 
very low compactness and large maximum Lorentz factor of 
the electrons have to be postulated, in order for GLAST to 
make a similar detection. 

These results are consistent with the hypothesis that 
there always exists a spectral component due to a non- 
thermal proton distribution in black hole binaries. This com- 
ponent peaks at 7-rays and can contribute between 0.5 to 
10% of the bolometric luminosity. In the ultra-soft state, this 
component is visible for E > 5 keV. During the soft and hard 
states, the emission is detectable only when E > 0.8 MeV, 
since other spectral components dominate at lower energies. 
This hypothesis can be verified using future observations by 
GLAST. 

The scheme does not include scattering of high energy 
photons with thermal electrons. Although the Klein-Nishina 
cross section decreases with photon energy, such scatterings 
can be important especially when the Thompson optical 
depth is significantly greater than unity, a case more per- 
taining to the soft state. In fact, the inability of the present 
model to explain the hard X-ray emission during the soft 
state, may be an artifact of this assumption. If that is true, 
then it may be possible that only protons are accelerated in 
black hole binaries and not electrons. This is theoretically 
appealing given the complexities of accelerating electrons to 
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high energies in a region where inverse Compton cooling is 
efficient. However, several sophisticated modification of the 
present scheme have to be undertaken before this specula- 
tion can be tested. 

Evidence of non-tlicrrnal protons in black hole binaries, 
would shed light on the energy dissipation process that occur 
in the inner regions of the accretion disk. These energetic 
protons could also be the origin of the outflows/jets that arc 
observed in many of these systems. Such insights may finally 
lead to a comprehensive physical picture of these enigmatic 
sources. 
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